Non-locality of non-Abelian anyons 
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Topological systems, such as fractional quantum Hall liquids, promise to successfully combat environmental 
decoherence while performing quantum computation. These highly correlated systems can support non-Abelian 
anyonic quasiparticles that can encode exotic entangled states. To reveal the non-local character of these en- 
coded states we demonstrate the violation of suitable Bell inequalities. We provide an explicit recipe for the 
preparation, manipulation and measurement of the desired correlations for a large class of topological models. 
This proposal gives an operational measure of non-locality for anyonic states and it opens up the possibility to 
violate the Bell inequalities in quantum Hall liquids or spin lattices. 



Quantum mechanics is a non-local theory: it allows for cor- 
relations between distant systems that cannot be explained 
in terms of a local preparation. Many believed that non- 
locality was due to incompleteness of quantum theory. Ein- 
stein, Podolsky and Rosen (EPR) in their seminal work [1] 
aimed to resolve this by introducing local hidden variables. 
Their values would complement the information supplied by 
quantum mechanics, thus restoring locality. Bell inequali- 
ties aim at validating or rejecting this view from experimental 
data [1]. To date, unlike local hidden variable (LHV) theories, 
the predictions of quantum mechanics have been consistent 
with all Bell tests. 

Particle exchange gives a striking example of non-locality 
in quantum mechanics. For bosons and fermions, one can deal 
with the exchange interactions by imposing non-local con- 
straints on the form of the wave function. More generally, 
the wave function can transform in a nontrivial representation 
of the fundamental group of configuration space when parti- 
cles are adiabatically exchanged [ ]. For planar systems, this 
group is the braid group and particles transforming in nontriv- 
ial braid group representations have been dubbed anyons [4]. 
Anyonic exchange interactions are topological in nature and 
do not change on variation of the distance between the parti- 
cles, or the metric of the spacetime manifold. One mechanism 
for such interactions is the Aharonov Bohm effect [5]. 

It is natural to ask if the non-local correlations of anyons 
can in principle be explained using local hidden variables. 
In this paper we answer this question for a number of anyon 
models, by constructing a Bell test for anyonic degrees of free- 
dom and showing violation of the associated Bell inequalities. 

Anyons can be split into two main categories, they can be 
Abelian or non-Abelian. Given labels {aj} for the different 
types of anyons we assign fusion rules that determine the out- 
come of bringing two anyons together, a; x aj — Y.k^a^aj'^k- 
Here A'^^ G N counts the number of ways of combining a and 
b to obtain c. Non-Abelian anyons have Y.c^ab — ^ ^'^^ some 
pair a,b, while in the Abelian case, the labels of the fused 
anyons determine a unique outcome which can be determined 
in a unique way. 

In a physical system with anyons, the low energy part 



of the Hilbert space can be thought of as a tensor product 
y{ = .?^iocai <8) .?^non-iocai, where the first factor describes lo- 
cal degrees of freedom, which we will ignore, and the sec- 
ond describes topological degrees of freedom associated with 
the anyons. These topological degrees of freedom may arise 
as a result of nontrivial topology of the space supporting the 
anyons. For Abelian anyons, this is in fact the only possibility; 
in the Abelian toric code models [d] for instance, the non-local 
degrees of freedom are described by elements of the first ho- 
mology groups of the surface with finite group coefficients. In 
principle, one probe non-local correlations in these topologi- 
cal degrees of freedom, but the observables involved would 
need to be non-local themselves [ /]. 

For non-Abelian anyons, even on a contractible surface, 
there are non-local degrees of freedom associated with the 
different fusion outcomes. A number of proposals have been 
made on how the associated quantum numbers, or topolog- 
ical charges, might be measured by interferometry, see [8], 
[y] for further references and [10] for an overview of the 
measurement theory. We will not go into the details of in- 
terferometric measurements here, but rather just assume that 
we can do projective measurements onto the various fusion 
channels. The non-local Hilbert state space of n anyons 
(fli,fl25 ■ ■ with total charge c has dim(y/non-iocai) = 

I/.. .b, b„_, N^>\cnKUKla, ■ ■ ■Nl_,a., ■ This Hilbert space 

usually does not admit a tensor product structure, e.g. the 
dimension could be prime, and thus does not obviously fit 
the usual paradigm for tests of non-locality. Nevertheless, 
we show that topological interactions can indeed be used to 
demonstrate non-IocaUty in the EPR sense. In order to do 
this, we consider two classes of anyonic theories: the SU{2)k 
models, including the Fibonacci model [9], and a model based 
on discrete gauge theory [11, 12]. These models are impor- 
tant both for their potential to process quantum information 
fault tolerantly and for their viability for experimental real- 
ization. For these cases the fusion spaces are at most one 
dimensional, i.e. A^^^^ < 2 for all {a,b,c). Non-commuting 
measurements project onto different ways of combining par- 
ticles a,b,c to yield d. Measurement bases are labelled by 
the intermediate products x and x' obtained by fusing a,b,c 
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and are related by the recoupling formula: \{ab)c d;x) = 

To date, most Bell tests have been performed on entan- 
gled light beams, but there is certainly an interest in showing 
that material media can be used to demonstrate non-locality. 
Experiments involving a photon and an atom, two atoms, or 
even kaons have been proposed or even carried out [13]. The 
schemes we are proposing contribute to this effort of using 
ever new media for Bell tests. As we shall see, some of them 
could be implemented in a fractional quantum Hall liquid [M], 
while others can be associated with arrays of Josephson junc- 
tions [ I ] or atoms in optical lattices [15]. 

In order to build intuition for the anyonic case we describe 
the general framework by employing distinguishable spin- 
1/2 particles with non-trivial fusion properties. The fusion 
rules are interpreted as the angular momentum decomposi- 
tion of tensor products of vector spaces. Consider a system 
divided into two spatially non-overlapping subsystems A and 
B, conveniently labeled as Alice and Bob, each one possess- 
ing three spin- 1/2 particles, as seen in Fig. 1. First, we per- 
form a joint measurement on the total spin 5tot = Lj ^nd 
post-select the Stot = outcome that has state space dimen- 
sion five. Second, we define a set of measurement operators 
V^iiJiiJis^'^leh where T,,, = {si + Sjf- - 1. The eigen- 
values of Y, j are +1 in the triplet space and —1 for the singlet 
and the operators y^f^) act on the subsystems A{B). The op- 
erator pair within A or B is non-commuting but [T^^- , T^^,] = 0. 
Consider the expectation value of the operator 
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(1) 



For a classical theory, even in the presence of local hidden 
variables (LHV) [16], the Bell inequahty for W is | (H')lhv | < 
2. This can be derived straightforwardly as follows [17]. 
Assume independence of the two subsystems (locality) so 
that the joint probabilities for pairs of outcomes is just the 
product of the individual probabilities which could depend 
on a hidden variable A,, drawn from a fixed distribution 
p(k). For the above quorum of observables with outcomes 
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Z2 3,'n45,'«f g} e ±1 we have 
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('«2,3+<2)<5-('«2,3- 

Hence, in the LHV model, the outcomes must satisfy 

iwiHvl = \idXpiX){iTUX)rl,iX) 

+TUX)T^,,{X)))\ 
= \JdXp{X){mj,{X)+mi,{X))mi,{X) 

-{n4^,{X)~mUX))ml,{X)\ 
< 2 

Quantum mechanically, the maximum value of | {W) \ is ob- 
tained for eigenstates of W with maximum eigenvalue, i.e. 

< Vl. For arbitrary operators in Eq. (1) that have the 
same commutation structure and square to 1, quantum me- 
chanics satisfies Tsirelson's inequality [IS], < \/8- Our 
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FIG. 1: A Bell type measurement on six particles. First a joint mea- 
surement (large oval) on all six particles is made and the result kept 
if the total charge (or spin) is zero. Alice performs measurements of 
total charge on pairs 1,2 and 2,3 and Bob performs measurements 
on pairs 4,5 and 5,6. For some quantum states the correlator (W) 
exceeds the bound set by local hidden variable theories. 



aim is to find a violation of the classical upper bound in the 
subspace of states with ^tot = 0. Note that our protocol allows 
for measuring correlations without the need of a shared refer- 
ence frame between Alice and Bob [ i >] thus giving a simple 
and unambiguous test of Bell inequalities. In the anyonic case 
treated below the operators T, j also have eigenvalue — 1 when 
the fusion outcome is the vacuum and +1 otherwise. 

There are two natural orthonormal bases for a three particle 
system based on the two different orders of fusing the three 
particles. These are graphically represented by fusion trees in 
Fig. 2b. The unitary transformation that describes the change 
from one of these bases to the other is given by the so called 
F matrices. The basis change for three particles with charges 
{a,b, c) fusing to d is given by a matrix denoted F^j'^ For the 
case of SU{2) these matrices just describe angular momen- 
tum recoupling and their matrix elements are the Wigner 6- j 
symbols. For six particles with total spin 0, we get four nat- 
ural product bases from the two pairs of bases for each triple. 
The contributing particle labels are the spin values {0, j , 1 , | } 
and the only relevant F-matrix for changing between bases is 

F^, I , which, for the SU{2) case is given by 
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in the basis given by fusion trees with intermediate spins and 
1 . We find the following set of orthonormal states for the six 
particle system, in the local fusion basis as defined in Fig. 2, 

|c^o) = |0')a|0)b, |^i) = |1')a|0)b, |(^2) = |0')a|1>b, 
|c^3> = |1'>a|1>b, |^4> = |l(ii))A|l(ii))i., 

(2) 

where \x) = \x{^, |)) and jjc') = T^x^^jI^)- ^P^^ components 
we have |(j)o) = |^")i,2 » |'P")3,4 ®V")5,6, where |^-) = 
(I Ti) ^ I iT))/V2, so the state |(|)o) has three adjacent singlet 
pairs. Notice that in order to have trivial total charge the local 
bases occur in pairs that share the same label p as defined in 
Fig. 2. 

In the basis { 1 0) |0) , 1 0) 1 1 ) , 1 1 ) |0) , 1 1 ) 1 1 ) } U |(^4) , we have 
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FIG. 2: The state space of anyons in our protocol represented as 
fusion trees, (a) An arbitrary state of six a type anyons witli trivial 
total charge expanded in terms of fusion outcomes local to A and 
B. In the models considered here, P is its own antiparticle, but it is 
straightforward to generalize, (b) A "local" fusion basis satisfying 
T^3|x(a,p))A = ±|x(a,P))A andXf 2|y(a,p))^ = ±|x'(a,p))A for 
X a vacuum state or a particle and similarly for B. For each subsystem 

the bases are related by an F move: |x'(a,P)) = Lv(^aaa) v'l'''(*^' P))- 



The Hilbert space splits into different sectors labeled by P, 
which are conserved by the action of W. There is a four di- 
mensional sector with (3=1/2 and a one dimensional sector 
with P = 3/2, containing |(|)4). No Bell violation can occur 
in the p = 3/2 sector, since the measurement operators com- 
mute in that sector Maximally Bell violating states are thus 
orthogonal to |(|)4). For the one parameter family of states 



|K«)> = ^(|^0) + |^3»- 
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(|c^i)-|^2», (4) 



with (^1 < fl < 1) we plot the expectation value (W) seen 
in Fig. 3. The maximal violation ((Vt')max.min = 



±2.6458) is obtained for fl± = (7±2\/7)/14. 

Consider now a two dimensional system with quasiparticle 
excitations described by SU{2)k Chern-Simons-Witten theo- 
ries. The corresponding fusion rules satisfy the addition of 
angular momentum with the constraints j\ x 72 j only if 
ji, ji ,jl^k/2 and j\ + 72 +j<k. It is quickly verified that for 
A: > 3, the total charge zero sector of six particles labeled by 
spin 1/2 again has five states, labeled by the same fusion trees 

as in the SU{2) case. The F matrices will differ, but for our 

_ 1 

purposes, the only relevant recoupling is still , , . Comput- 

III 

ing the quantum 6- j symbols, we find (see for instance [20]) 
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where the quantum integers are defined as [m 
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for m integer. For the SU{2)k theories, q = e^. In the limit 
k ^ then [m]c/ m. As before, we can label the states in 
the local fusion basis, as in Eq. 2, but with the appropriate F 



matrix. The state |(|)o) is obtained by creating spin-1 /2 parti- 
cle anti-particle pairs at positions (1,2), (3,4), (5,6) out of the 
vacuum. For the one parameter family of states \r{a)) we find 
maximal violation at 
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and at fl_ 



1 — a-., where 
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It is possible to verify that k ^ °o corresponds to SU (2). A 
qualitative difference between anyonic systems and the spin 
systems discussed before is that, while the z-components of 
the spins of all particles are in principle measurable, there 
are not necessarily any observables associated with the z- 
components of the 'q-spins' of the anyons. Only SU{2)g 
invariant quantities, such as the total q-spins of groups of 
anyons, can be observables, or at any rate topologically pro- 
tected observables. This can be traced back to the superse- 
lection rule that says that the total q-spin of all anyons to- 
gether must be trivial. If it were possible to measure the z- 
components of every anyons' q-spin, then the state obtained 
would no longer be invariant under SU{2)q. In fact, a similar 
rule would hold for confined particles in gauge theory and so 
for a better analogy, one may think of the SU (2)^ invariance 
as begin closer to an SU (2) gauge symmetry rather than spin. 

It was shown by Freedman et al. [21] that the anyonic the- 
ories with ^ > 3, A: 4 are universal for quantum computa- 
tion. Hence, for those theories, the Bell violating states can be 
obtained by topological braiding operations alone acting, for 
example, on the fiducial state |(|)o). We now check to see if it 
is possible to generate a state which violates the inequality for 
the k^2 case. The SU (2)2 anyons are believed to exist in the 
V = 5/2 plateau of the fractional quantum Hall effect [? ] up 
to charge factors that affect the Abelian part of braiding. They 
come in three varieties, the vacuum, 1, the fermion, \|/, and 
the non-Abelian anyon, a, that satisfy the non-trivial fusion 
rules, 0x0=1+ \|/, ax\|/ = a and \|/ x \|/ = 1 . The counter- 
clockwise exchange of two o particles, which fuse to either 1 
or \|/, results in the matrix evolution R= 1 ® ; expressed in the 
basis labeled by the fusion channels The state evolu- 

tion produced by the exchange of particles with no immediate 
fusion channel is found by employing the recoupling matrix 
F. Expressed in the basis {|x(o,o))^|3'(o,o))B;A:,y G {1,\|/}}, 
we have the following representation of the generators of the 
braid group 
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(6) 
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where Bj results from the exchange of j and j +1 particles in 
a counterclockwise manner As a simple initial state we can 
consider |(|)o) that is produced from (1,2), (3,4) and (5,6) 
pairs created from the vacuum. 
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FIG. 3: The expectation value of the Bell witness ly as a function 
of the amplitude of mixing for the total charge zero states \r{a)) in 
Eqs. (4). The yellow, blue, and red lines correspond to an SU{2), 
50(3)3, SU (2)2 theory with six spin-1/2, a, x particles, respectively. 
The shaded region corresponds to states which violate the inequality 
derived for local hidden variable models. 

The braid group generators, Bj, are in the Clifford group, 
so we cannot generate a dense set in SU (4) by braiding 
alone[31]. But can we still obtain Bell violating states? In 
Refs. [24, 25] a LHV model was introduced for a pair of 
qubits that exactly reproduces the set of allowed operations 
in the present model. There are two distinct configurations of 
shared vacuum pairs (up to relabeling of particles by Alice or 
Bob) both of which can be obtained from |(|)o) by braiding. 
Hence it is not possible to build Bell violating states starting 
out from three shared vacuum pairs using topologically pro- 
tected operations alone. 

Despite the impossibility of producing a Bell violating state 
from |(|)o) by topological gates, one can in fact straightfor- 
wardly obtain a maximally Bell violating state using non- 
topological gates[32]. Let us employ the non-Clifford gate 
D = e^'^^' (g) 1. This can be implemented by bringing the 
2 and 3 O anyons nearby, thus shifting the energy of the 
fermionic fusion channel [27] such that a relative phase 
gm/4 jg accumulated on that channel. From these opera- 
tions one can build the controlled phase gate in the follow- 
ing way CP = e™/'*B2BiB2B3 '^2 ^ simple search- 
ing algorithm provides us with the sequence that produces 
|r(fl_)) = -CPBj.B4DB2Bj,\<:^o), with (W) = -2\/2, thus sat- 
urating the Tsirleson bound. 

In fact, it is indeed possible to realize a maximally violat- 
ing state without braiding at all. Consider a state given 
by a distribution of singlet pairs on (1,6), (2,5) and (3,4). 
This state |(|)q) is related to the fiducial distribution of pairs by 
the following braid word |^'o) - ^(|0')a|0)b + |1')a|1)s) = 

B^^ B^^ BsB/iBT,B2\^Q) ■ This state would be maximally vio- 
lating if we could measure in arbitrary local bases. For our 
fixed measurement quorum |(|)q) is related to the maximally vi- 
olating state by e^'^°' (g) 1|(|)q) = \r{a+)) and hence it suffices 



to implement the local unitary e^'v°^ = e^'f ^'g^'j'^'e'f ''"^ on 
Alice's side to obtain a Bell violation. The z rotation is sim- 
ply achieved by bringing anyons 2 and 3 near each other as 
above, and similarly the x rotation is performed by pushing 1 
and 2 together Note that the maximal Bell violation in these 
two constructions actually saturates the Tsirelson inequality, 
making the SU (2)2 case at the same time 'maximally quan- 
tum mechanical' and 'topologically classical'. 

Let us turn now to Fibonacci anyons from the 5(9(3)3 the- 
ory. This is the theory obtained from SU (2)3 but using only 
integer spin particles: the vacuum 1 and the non-Abelian 
anyon 1, with non-trivial fusion rule T x T = 1 + T. All par- 
ticles are their own anti-particles and the quantum dimensions 
are Ji = 1 and (ix = (j) = ( 1 + \f5)/2. The relevant recoupling 
matrix is 

expressed in the basis of 1 and T. The dimension of the 
topological Hilbert space of m + 1 type T anyons with total 
charge zero is /,„, the mth Fibonacci number, hence there are 
five states in the fusion space. These states can be decom- 
posed into superpositions of products of local basis states as 
in Eq. 2 where |0) = |1(t,t)) and |1) ^ |t(t,t)) and {(^4) ^ 
|t(t, 1))/i|t(t, 1))b- The state |(|)o) is the state obtained by cre- 
ating type T particle anti-particle pairs on (1,2), (3,4), (5,6) 
out of the vacuum. For the one parameter family of states 
\r{a)) we find the same maximal violation as in Eq. 5 for 

SU{2)3: (W) = ±2v/-7+4v^« ±2.7887. This is not very 
surprising, considering that the SU{2)3 theory is equivalent 
to the product of the Fibonacci theory and an Abelian theory 
with Z2 fusion rules (see for instance [28]). The action un- 
der braiding is represented by the matrix /?xt = e'^''/^ g'7it/5 
expressed in the basis {1,t}. We obtain the following repre- 
sentation of the generators of the braid group expressed in 
the basis {|0)|0),|0)|1),|1)|0),|1)|1)}U 1^4): 

Bi = [F/?xx-F-'®l2]®(e'^"/^) 
B2 = [^xT<^l2]®(e''''/^) 

Bi = O' [e'4V5 © V5 ^ e jjji ^8 ^ ] O (7) 

B4 = [l2(gF7?xTi^-^]©(e"'^/5) 
B5 = [l2®^Tx]©(e'^''/^) 

where O maps the product basis to the basis {|(|);)}j^0' ^^'^ 
M = (Fxxx)"'/?xT^'xxT- A length 25 braid word produces a Bell 
violation: = [B3B4 '^r'^a '^2 'l^l^l^o) with (1'|W|1') = 
2.5310. 

In the models above, measurements by Alice and Bob had 
two outcomes for the two fusion products of the anyons. To 
accommodate more outcomes we can use higher dimensional 
Bell witnesses [2 ']. We demonstrate how this works in an- 
other anyonic model with excitations in one to one corre- 
spondence with irreducible representations of a Hopf algebra, 
D{G), the quantum double of a finite group G [1 1, 12]. The 



5 



particles can carry electric and magnetic charge and are la- 

beled by n^(^| where [a] denotes a conjugacy class of G 

which labels the magnetic charge, and /?(A^[oi]) denotes a uni- 
tary irreducible representation R of the centralizer of an el- 
ement in the conjugacy class [a], which labels the electric 
charge. The dimension of the carrier space for each irre- 
ducible representation, which equals the quantum dimension 
la] 

of the particle n^A, ^ i& d, 



«(A'[al) 



|[a]||7?(A?[„])|. The quan- 



tum dimensions satisfy the sum rule ^('^^"Ivf |)^ ~ l^l"' 
focus on the simplest non-Abelian finite group, 53, the group 
of permutations on three objects. Elements of S3 are orga- 
nized into three conjugacy classes: [e] ~ {e} the identity el- 
ement, [t] = {tQ,ti,t2} the transpositions, and [c] = {c+,c_} 
the cyclic permutations. The 8 irreducible representations for 
0(83) are 



n 



vacuum) 

2, 3 (pure magnetic charges) 
1 , 2 (pure electric charges) 
Hpl, 11^^,11^1 d = 2,2,3 (dyonic combinations) 



' n" d 



(8) 



A complete derivation of the fusion rules for this model is 
given in [30]. In the toric code realization of these anyon mod- 

els, the quantum dimensions d^^^^ ^ = |[a]||/?(A^[oi])l actually 
count local degrees of freedom associated with the anyon. In 
the discrete gauge theory context, these degrees of freedom 
are also present in the description of the system, but some of 
them are gauge. A single particle's electric charge and mag- 
netic charge can always be measured locally (or at least within 
a region of size characteristic of the particles), by braiding 
with other locally prepared charge pairs and measuring the 
outcome of fusion of the pairs. Truly non-local properties are 
contained in the fusion space. To explore this we pick a fu- 
sion subalgebra of 0(83) : {n^l,n^^L,n^^j} which we label 
for convenience {l,A,<t>}. The non trivial fusion rules are 

AxA=l, Ax4> = 4>, X <I> = 1 +A + 4>. 

These fusion rules are the same as the fusion rules for the 
representations of 53 itself and also the same as the fusion 
rules of the integer spin sectors of 8U{2)4. The particles are 
their own anti-particles. The magnetic charge <!> with quan- 
tum dimension 2 carries non-Abelian statistics and the fusion 
of « such particles gives: 4>^" = 1(2"-' + (-!)")(!+ A) + 
1(2" + {-1)"-^)<P. As before, we will work in the superse- 
lection sector with total trivial charge. The smallest number 
of particles in this sector that could hope to violate a Bell in- 
equality should have fusion space dimension > 4. If we are 
to pick measurement operators for Alice and Bob that mea- 
sure total charge on pairs of <I> particles and we want two non 
commuting operators on each side then we require at least six 
particles in total. Exactly six particles suffices, giving Hilbert 
space dimension eleven for the vacuum sector. 



Either by using the representation theory of 0(83), or by 
solving the pentagon and hexagon equations directly, we find 
the following recoupUng and braid matrices, expressed in the 
basis {1,A,4>}, 



F = F^c^^ — 
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\;R = R^<s> = 1 
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We notice immediately that R has eigenvalues ±1, so that we 
will end up with a representation of the permutation group 
when 'braiding' the anyons. Nevertheless, these anyons are 
not bosons or fermions, since this representation is non- 
Abelian. The fact that we have a permutation group repre- 
sentation does signal the fact that braiding in this theory is not 
universal for quantum computation. This is in fact a general 
property of braiding in discrete gauge theories. 

A basis of the eleven dimensional vacuum sector of the six 
anyon Hilbert space can be given in terms of superpositions 
of products of local basis states, as defined in figure 2, 

{L.F'^\y{<P,<P))A\l{<P,<P))BXyFl\y{'P,<P))A\l{<P,<P))B 
LrFl\yi<f,<f))A\l{<f,<f))B,LyF^'\y{<P,<P))A\A{<f,<P))B 
I, Fl\yi<P, <P))a \A{<P, <P))B,Ly Fl\yi<P, <P))a \A{<P, <P))b 
LyFl'\yi<P,<f))AM<f,<f))B,i:,Fl\yi<f,<P))A\<P{<P,<P))B 
^,Fl\yi<P,<f))Am<P,<f))B,M'f,A))A\<f{<P,A))B 

\<p{<p,i))AM<p,i))B} = {m]%- 

The state |(|)o) is obtained by creating type «t> particle anti- 
particle pairs on (1,2), (3,4), (5,6) out of the vacuum. Each 
such vacuum magnetic charge pair is written: = 
75 ( |c+ , c_ ; (/, 7)) + |c- , c+; (/, j))). 

Now in analogy to the cases studied for 8U{2)k, we could 
look for a Bell like inequality but using measurement oper- 
ators with three outcomes. Let Alice have one operator T'^ 2 
which measures the outcome of total charge for particles 1 
and 2 with outcomes {l,A,<t>} and another, non commut- 
ing operator, T23 that measures total charge for particles 2 
and 3 with outcomes {1,A,<I>}. In other words, is a 

measurement in the basis {(faaa)^b(0C5p))A} with outcome 
'"f 2 ~ y ^ {^7A,<P} and Tj 3 is a measurement in the basis 
{|v(a, P))a} with outcome m2 3 =y. Similarly, let Bob have 
two measurement operators, T45 that measures in the basis 

{(^'laa)'''b(oc,P))B} with outcome 5 = y, and Tf g which 
measures onto the basis {|}'(oc, P))b} with outcome = y. 
Now {Fl^^)l = 5x,<i>5y,<i, = {F^^^)x, so in the subspace of 
{\i^g), |*j)io)}, the measurement operators all commute. These 
states cannot yield a Bell violation and we can focus on the 
states in the 9 dimensional orthogonal subspace which is iso- 
morphic to the Hilbert space of two three dimensional parti- 
cles (qutrits). 

In Ref. [29] it was shown how to construct Bell inequalities 
for bipartite systems of equal but arbitrary finite dimension. 
In particular for two qutrits the authors introduce the witness 
I3 which for all LHV theories satisfies |(/3) | < 2, whereas for 
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quantum mechanical systems \{h)\ < 4. To simplify notation 
let us introduce the projectors Ky = |3'(<I>,<I>))(3'(<t>,<t>)| and 
%y = F^y{(i>,<i>)) {y{ci>,ct>)\F . For the quorum of observables 
above, the Bell witness It, is; 

+7tA '8' tia + ft* ^ — 1 <8 Sa + Sa (8) ft* + ft(I> ® ft 1 

— Tti (8) fti + 7tA g) ftA + <8l ft* — Jtl (g JtA + 'Ia 71$ 
+7t<j> (g) Tti — ftA CS" Ttl + ft* ® TIa + ftl ® Tt* 

+2|4>(4>, A)) (4>(4>, A)| ® |4>(4>, A)) (4>(<I>, A) | 
+2|4>(4>, 1)) (4>(4>, 1)1 ® |4>(4>, 1))(4>(4>, 1)1 

(9) 

The state with the largest violation has (/s) = —2.5216. 

We obtain the following representation of the generators 
for -Be expressed in the basis {|x(<I>,<I>))a|3'(<I>,<I>))b;x,v e 
{1,A,4>}}U|4>(4>,A))a|4>(4>,A))bU|4>(4>,1))a|4>(4>,1))b}: 



Mi 



K 



Mf 



at least six anyons shared between two parties arises because 
each party needs three anyons in order to have two non- 
commuting topologically protected observables. It is possible 
this could be reduced using a shared resource which fixes a 
common gauge, akin to using a shared reference frames to re- 
veal non-locality in mode entanglement with bosons [19]. The 
size of the maximum violation depends on the recoupling ma- 
trices F and the ability to generate Bell violating states begin- 
ning from 3 vacuum charge pairs depends on the power of the 
braiding operations. It is intriguing to ask whether one could 
find intermediate anyonic theories which have the power to 
generate Bell violating states by topologically protected gates, 
but are not universal for topological quantum computation. 
Finally, it would be very interesting to have an experimen- 
tal demonstration of (some of) the schemes presented here. 
We believe that the required effort would not be significantly 
higher than that necessary to perform non-abelian interferom- 
etry. 
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B4 = [h®FRF-'](SR%^(SR%^,B5^[h®R] 



where O maps the product basis to the basis {|(|);)}]|io, 
M = F^^RF. Here = InVy, so we have the permutation 
group 56, as mentioned before. We compute the action on 
states consisting of vacuum magnetic charge pairs. The state 
|(|)o) = |<I>,<I>;(1,2))|4>,<I>;(3,4))|4>,4>;(5,6)) is the fiducial 
state, and the other distinct configuration of vacuum magnetic 
charge pairs is |4>,4>; (1,4)) |<I>,4>; (2,5)) |4>,4>; (3,6)) = 
l^l^o), hence it suffices to consider the orbit of |(|)o). 
An exhaustive search through 6! = 720 braid words corre- 
sponding to all distinct permutations in Se finds that, while 
(73) is not constant under braiding, we do find that in all cases 
1(73) I < 2. Hence we require some operation beyond braiding 
to produce a violation of LHV under our protocol. Fven if 
we restrict to non-topologically protected operations that just 
involve interacting pairs of particles, we can indeed produce 
a Bell violating state. Consider the family of states |(|)') = 
£>3,4(oci , a2)£>i,2 (0C3 , a4)7)2,3 (as , ae)BiBiBo,B2B^B4 |(j)o) 
where D, j(a, p) is the non-topologically protected gate 
obtained by bringing anyons / and j of type <1> nearby 
each other and allowing them to interact for a time such 
that the fusion channel <I> x <I> A accumulates a phase 
e'" and the fusion channel <I> x <J> ^ <I> accumulates a 
phase e'P. Optimizing |(73)| over the interaction phases, 
we find a violation ((|)'|73|(|)') = 2.0512 for the angles: 
ai = 0.7943, a2 = 0.3989, a3 = 3.5531, a4 = 0.9257, as = 
— 0.8525, ag = 0.1036. No systematic attempt was made to 
optimize the violation over other braid words and it likely 
stronger violations could be found. 

We have described a protocol to reveal non-locality in sev- 
eral classes of non-Abelian anyonic theories. The need for 
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